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Abstract 

Non-Markovian transport equations for nuclear large amplitude motion are 
derived from the collisional kinetic equation. The memory effects are caused 
by the Fermi surface distortions and depend on the relaxation time. It is 
shown that the nuclear collective motion and the nuclear fission are influenced 
strongly by the memory effects at the relaxation time r > 5 • 10~ 23 s. In 
particular, the descent of the nucleus from the fission barrier is accompanied 
by characteristic shape oscillations. The eigenfrequency and the damping of 
the shape oscillations depend on the contribution of the memory integral in the 
equations of motion. The shape oscillations disappear at the short relaxation 
time regime at r — ► 0, which corresponds to the usual Markovian motion in 
the presence of friction forces. We show that the elastic forces produced by 
the memory integral lead to a significant delay for the descent of the nucleus 
from the barrier. Numerical calculations for the nucleus 236 U shows that due 
to the memory effect the saddle-to-scission time grows by a factor of about 
3 with respect to the corresponding saddle-to-scission time obtained in liquid 
drop model calculations with friction forces. 

PACS number: 21.60.Ev, 25.85.Ca 



1 



Typeset using REVlp^X 



I. INTRODUCTION 



The dynamics of a nucleus undergoing fission can be studied in terms of only a few 
collective variables like nuclear shape parameters Such kind of approach is usually as- 
sociated with the liquid drop model (LDM) and its extensions and is acceptable for a slow 
collective motion where the fast intrinsic degrees of freedom exert forces on the collective 
variables leading to a Markovian transport equation. An essential assumption is that the 
LDM provides a good approximation for a smooth part, E pot , of the collective potential 
energy, E pot , and can be then used for the quantum calculations of E pot within Strutinsky's 
shell correction method obtaining E pot = E pot + 5U, where SU is the shell correction. 
On the other hand, it is well known that the LDM is not able to describe some strongly 
collective nuclear excitations such as the isoscalar giant multipole resonances. It is because 
the LDM ignores the important features of the nucleus as a Fermi liquid. The collective 
motion of the nuclear Fermi liquid is accompanied by the dynamical distortion of the Fermi- 
surface and the smooth energy E pot is subsidized by an additional contribution, -Epo^F, 
which is caused by the dynamic Fermi-surface distortion effect and is absent in the stan- 
dard LDM HIJ. We point out that the energy -Epot.F is a smooth quantity (in the sense 
of the shell correction method) and it can not be recovered by taking into consideration 
the quantum shell corrections to the adiabatic (static) potential energy deformation. This 
situation becomes more clear in the limit of the infinite Fermi system with A — > oo, where A 
is the number of particles. The shell correction SU disappears at A — > oo. In this case, the 
adiabatic collective energy E pot = E pot , caused by a change of the particle density p with 



respect to its equilibrium value p eq , determines the first sound velocity C\ = ^jK/9m, where 
K is the incompressibility coefficient given by 



d 2 (E pot /A) d 2 (E pot /A) 2 



However, it is well known |||7j that, in a cold Fermi liquid, the first sound velocity C\ is 
relevant only in the limit of strong interaction, i.e., at \Fq\ ^> 1, where Fq is the Landau 
parameter in the quasiparticle scattering amplitude. For the normal nuclear matter we 
have rather the value of Fq ~ 0, which can be derived from the Skyrme forces ||. The 
corresponding sound velocity (zero sound velocity) c exceeds the velocity c\ by a factor of 
about \f?> at F ~ [[7|. This difficulty is overcome if the additional contribution E pot:F to 
the potential energy E pot is taken into account |J. Thus, the smooth energy E pot ^ F , caused 
by the Fermi surface distortion effect, is a necessary ingredient of the dynamics of the nuclear 
Fermi liquid. It is absent in the adiabatic deformation energy E pot derived by the traditional 
LDM. 

The equations of motion for the nuclear Fermi liquid can be derived from the collisional 
kinetic equation ||. In general, the corresponding equations of motion are non-Markovian 
T0| , |lT]| . The memory effects appear here due to the Fermi-surface distortion and depend on 



the relaxation time |12| , p!3| . The Markovian dynamic is achieved in two limiting cases only: 
(i) short relaxation time limit which corresponds to the first sound propagation in infinite 
Fermi liquid. In fact, this limit is realized by the nuclear LDM, (ii) infinite relaxation time 
limit which corresponds to the zero sound propagation with a strong renormalization of the 
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sound velocity and the deformation energy with respect to the ones in the LDM. The non- 
Markovian-Langevin equations of motion for macroscopic collective variables were earlier 
derived in [H] and used for the small amplitude dynamics p|,|T5f and for some aspects of 
the induced nuclear fission and the fission rate problem [|l^ . 



The main purpose of the present paper is to apply the non-Markovian dynamics to the 
descent of the nucleus from the fission barrier. Starting from the collisional Landau- Vlasov 
kinetic equation, we suggest a new proof of the non-Markovian equations of motion for the 
nuclear shape variables which establishes a direct connection between the memory effects 
and with the dynamic distortion of the Fermi surface. In contrast to Ref. we do not 



take into consideration the random forces and only concentrate on the formation of both 
the conservative and the friction forces behind the saddle point to clarify the effects of the 
memory integral. In this aspect, our approach represents an extension of the traditional 
LDM theory of the nuclear fission [|H,|HHl9H to the case of the Fermi liquid and takes into 
account the important features of the dynamic Fermi surface distortion which are ignored 
in the LDM. 

The plan of the paper is as follows. In Sec. II we obtain the Euler-like equation of 
motion for the displacement field. This equation contains the memory dependent pressure 
tensor. Assuming that the nucleus is an incompressible and irrotational fluid and using 
the boundary conditions for the velocity field, we reduce the local Euler-like equation to 
the non-Markovian equations of motion for the shape variables. The transport coefficients 
and the memory kernel are derived through the solution to the Neumann problem for the 
potential of the velocity field. In Sec. Ill we study the dependence of the memory effects on 
the relaxation time for both the small amplitude motion near the saddle point and for the 
descent of the nucleus from the barrier to the scission point. Summary and conclusions are 
given in Sec. IV. 



II. NON-MARKOVIAN DYNAMICS OF NUCLEAR FERMI LIQUID DROP 

To derive the equation of motion for the shape variables, we will start from the collisional 
kinetic equation for the phase-space distribution function / = f(r, p;t) in the following 
general form 

TrJ + --Vrf-V T U-V p f = I[fl (1) 

at m 

where U = U(r,p;t) is the selfconsistent mean field and /[/] is the collision integral. The 
momentum distribution is distorted during the time evolution of the system and the distri- 
bution function takes the form 

/(r, p;t) = / sph (r, p;t) + ]T <J/,(r, p;t), (2) 

i>i 

where f sp h{ r ,P]t) describes the spherical distribution in momentum space and I is the 
multipolarity of the Fermi-surface distortion. We point out that the time dependent 
Thomas-Fermi (TDTF) approximation and the corresponding nuclear LDM are obtained 
from Eq. ([!]) if one takes the distribution function /(r, p;£) in the following restricted form 



3 



/tf(i", p;t) = / sp h(r, P',t) + 5fi=i(r, p;t) instead of Eq. (|), see Ref. p0[] . Below we will extend 
the TDTF approximation taking into account the dynamic Fermi surface distortion up to 
multipolarity I = 2 P,|5|,|2T|j . We will also assume that the collective motion is accompanied 



by a small deviation of the momentum distribution from the spherical symmetry, i.e., even 
in the case of large amplitude motion the main contribution to the distribution function 
/(r,p;t) is given the Thomas-Fermi term f^p(r,p;t) and the additional term 5fi=2(r,p;t) 
provides only small corrections. The lowest orders I = and 1 (which are not necessary 
small) of the Fermi-surface distortion do not contribute to the collision integral because of 
the conservation laws and the linearized collision integral with respect to small pertur- 
bation Sfi =2 (r, p;t), is given by 

/[/] = , (3) 

T 

where r is the relaxation time. 

Evaluating the first three moments of Eq. (JJ) in p-space, we can derive a closed set 
of equations for the following moments of the distribution function, namely, local particle 
density p, velocity field u u and pressure tensor P Ufl , in the form (for details, see Refs. [|20l . [2~2"| ) 

d 

T^P = -Vv(pUv), (4) 



—p' + -p_!^A = --P' (Q) 
dt v » + ' dt ^ r K ) 



where V = V(r,t) is the isotropic part of the pressure tensor 

P' Vil = P^(r,t) is the deviation of the pressure tensor from its isotropic part, V(r,t), due to 
the Fermi surface distortion 

Kfi( r J) = ~J (2ifh)^ Pv ~~ mUu ^ p » ~ mu v) 6 fi=2(r, p;t), (8) 

e pot is the potential energy density related to the selfconsistent mean field U as U = Se pot /Sp. 
The tensor A Ufl in Eq. (||) is given by 

2 

A VAI = V^Xm + V M x^ - 2^m V aXa, (9) 

where Xv = Xv( v t) is the displacement field related to the velocity field as u v = u u (r,t) = 
dxv(r,t)/dt. From Eq. @ we find the pressure tensor P^(r,t) in the following form 

P^(v,t) = P^(r,t ) - jf*df exp(^) ^Au,(r,t'). (10) 
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The tensor P'(r,to) is determined by the initial conditions. In the case of the quadrupole 
distortion of the Fermi surface, the tensor P' (r,t ) is derived by the initial displacement 
field Xv 

Assuming that the nucleus is an incompressible and irrotational fluid with a sharp surface 
in r-space, we will reduce the local equation of motion (|5|) to the equations for the variables 

q = {qi,q2, Qn} that specify the shape of the nucleus. The continuity equation ([|) has 

to be complemented by the boundary condition on the moving nuclear surface S. Below 
we will assume that the axially symmetric shape of the nucleus is defined by rotation of 
the profile function p = Y(z, {%(*)}) around the z-axis in the cylindrical co-ordinates p, z, <p 
P^pSfl . The velocity of the nuclear surface is then given by P5| 

N 

«s = X)^> ( n ) 

i=l 

where 



^ = (dY/d qi )/A, A = y/l + (dY/dzf. (12) 
The potential of the velocity field takes the form 



N 



(13) 



i=i 



where the potential field = 0j(r,g) is determined by the equations of the following Neu- 
mann problem |25| 



1 BY 

V 2 l= O, (nV0,) 5 = x ^-, (14) 

where n is the unit vector which is normal to the nuclear surface. 

Using Eqs. @ and fllCf ) with u v = V„0, multiplying Eq. (|) by V M 0j and integrating 
over r, one obtains 

EPvm + E ^ + fdt'eM^hait,*) q 3 (01 = ( i 5) 

j = i fc =1 vqk J to I uqi 

Here B>ij(q) is the inertia tensor 

Bij(q) = mpo j dstii fy, (16) 

where po is the nuclear bulk density. The adiabatic collective potential energy, E pot (q), does 
not contain the contribution from the Fermi-surface distortion effect and is given by 

E P ot(q) = J dr (e kin (r, q) + e pot (r, q)), (17) 

where ey m (r, q) is the kinetic energy of the internal motion of nucleons. The memory kernel 
Kij(t,t') in Eq. (|15|) is given by 
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n l3 (t,t') = 2 Jdv V(v,q(t')) (V„V„ ^ (r,g(t))) (V„V„ ^ (r,g(f))). (18) 

In Eq. fll5[) , we have omitted the contribution from the initial distortion of the Fermi 
surface caused by the pressure tensor P' (r,t )- The contribution from P' (r,i ) reflects the 
fact that the initial displacement field x^(r,t ) is switched on suddenly at t — t . The 
adiabatic force — dE pot (q) / dq in Eq. (|l5|) obtains then the additional contribution at t — to 
due to the initial distortion of the Fermi surface. However the corresponding force is absent 
if the initial displacement field ^(r,^) is obtained as a result of the previous evolution of 
the system at t < to- Below we will consider the descent of the nucleus from the fission 
barrier, i.e. assume the presence of the evolution of the system at t < to, and omit the 
contribution from the pressure tensor P^(r,i )- 

The displacement field Xv( r il) an d the potential field § { = (f>i(r,q) are determined by 
a solution to the Neumann problem (|TJ]). The displacement field can be also ob- 

tained using the Werner- Wheeler method f24|. In the cylindrical co-ordinates p,z,tp, the 
components of velocity field u z and u p in z and p directions are then approximated as 24 



u z = ^A l (z,q)q i , u p = P - ^Bi{z,q)qi. (19) 
i * \ z i 1) i 

The two unknown coefficients Ai(z, q) and Bi(z,q) are related to each other by means of 
continuity equation as 

B l (z,q) = - 1 -Y(z,q)^^l. (20) 

Requiring then that the normal velocity of the fluid on the surface should coincide with the 
normal velocity of the surface one can express the coefficient Ai(z, q) in terms of the profile 
function Y(z, q) as 

A l (z,q) = Y- 2 (z,q)—j z dz'Y 2 (z>,q). (21) 

We point out that in the case of irrotational flow the Werner- Wheeler method leads to a 
velocity field potential of quadrupole type 125 



<P(r,q) = j- q (2z 2 -x 2 -y 2 ). (22) 

A spheroidal figure presents the simplest example which is consistent with the velocity field 
potential (p2|). In this case q = q{t) is the elongation of the figure in units of the radius 
Ro = tqA 1 ^ of the nucleus and the equation of motion (|T5|) takes the following form 

m i i'= _ /Vexp(^MM') i (f). (23) 

oq oq Jt t 

Here, the mass parameter B(q) and the memory kernel K(t,t') are given by 

aw = ^^('4) and kM = WWV (24) 

where Kq = (4/5 m) tt p p 2 F -Rq and pp is the Fermi momentum. 
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III. NUMERICAL CALCULATIONS AND DISCUSSION 



Let us start from the one- dimension case and apply Eq. ( Ba ) to the large amplitude 
motion from the barrier point B to the "scission" point C in Fig. 1. Following the 
Kramers model |f26| , we will approximate the potential energy E pot (q) by an upright os- 
cillator (1/2)Cx£>m(<? — <?o) 2 with qo = 1 and an inverted oscillator Ef — (1/2)Cldm{q — <7/) 2 
which are joined smoothly as shown in Fig. 1 (see also Ref. f2?J). Let us consider, first of 
all, a small amplitude change, Aq, of the shape variable, q, near both the ground state at 
<Z ~ % — 1 with Aq = q — qo and at the saddle point at q ~ with Aq = q — qf- Linearizing 
Eq. (P5|), we will rewrite it as 

B ^-Aq = - k Aq - k /Vexp(— ) 4 A <?(0, (25) 



dt 2 ' Jt r dt 

where B = B = B(q =1), k — Cldm, k = kq if q ~ q and B = Bj = B{q = qf), k = 
—Cldm, K = Kf = Ko/q} if ? ~ qj- Differentiating Eq. (|25|) over time, we will look for the 
solution to Eq. ( p5|) in the form 

3 

Aq = Y / C i exp(X i t). (26) 

i=l 

Here the coefficients Cj are derived by the initial conditions. The eigenvalues Aj are obtained 
as a solution to the following secular equation 

(A 2 + i)(A + i) + |A = 0. (27) 

In the case of the zero relaxation time limit, r — > 0, one obtains from Eq. ( P?D a non- 
damped motion with A = ±\Jk/B, i.e., the time evolution is derived by the LDM stiffness 
coefficients Cldm or Cldm- In the opposite case of rare collisions, r — >• oo, the solution to 



Eq. (p7|) leads to a non-damped motion with A = ±y (A; + k)/B. In contrast to the previous 
case, the additional contribution, k, appears at the stiffness coefficient k + k because of the 
Fermi surface distortion effect. In the case of the nuclear Fermi liquid one has k 3> \k\ 0,0. 
This fact is important for the description of the nuclear isoscalar giant resonances f|J|. 
Considering a motion near the ground state at q ~ q = 1 with k = Cldm, one obtains 
from Eq. ( p7|) the quadrupole eigenvibrations with the eigenenergy 



hu 2+ = hJ ULD " + Ko « » J-^ « 64.5 MeV, (28) 




where e F = (9vr) 2/3 ^ 2 / 8 

mr\ = 34.73 MeV is the Fermi energy and we adopt r$ = 1.18 
fm. The result (|28|) coincides with the analogous one obtained earlier by Nix and Sierk Q 
and agrees with the experimental value of the energy of the isoscalar quadrupole resonance 
hull p ss 63 • A- 1 / 3 MeV. 

As can be seen from Eq. (|27D , the motion is damped for the non-zero and finite relaxation 
time r. In the case of small amplitude motion near the ground state, q ~ go, the solution to 
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Eq. ( f25|) at t — to ^ r takes the form of eigenvibrations with Aq(t) ~ exp(za;t), where the 
eigenfrequency u is derived by || 

^ 2 5 = Cldm + C'{u) - (29) 

where i?o = B(q = 1), the additional stiffness coefficient, C'(u), appears due to the Fermi 
surface distortion effect 

C'(u) = K kn(-^-) (30) 

VI — HOT ) 

and the friction coefficient 7(0;) is given by 

j(u) = k Re (-^—) = (4/m) vr R 3 Vo Re ( — L- ) . (31) 
\\ — iujt / \1 — iutJ 

Here, 770 = (1/5) po p 2 F r, is the classical viscosity coefficient @. 

Let us consider now the small amplitude motion (starting path for the development of 
the instability) near the saddle point, q ~ g/, at finite relaxation time. We have evaluated 
numerically the value of Aq from Eq. (|23| ) using the secular equation ( p7| ) and the initial 

conditions Ag(t ) = 0, A? (t ) = v and A g (t ) = 0. In Fig. 2 we show the result for two 
values of the relaxation time r = 3 ■ 10 _23 s and r = 4 ■ 10 _22 s. We have used the following 
parameters qj = 1.6, A = 236 and Tiujf = ?i\J\Cldm\/ Bf = 1.16 MeV. The initial velocity 
Vo was derived using the initial kinetic energy -Ekin.o = 0-/2)Bf v o — 1 MeV. In the case of 
the very short relaxation time, r = 3 • 10~ 23 s, the memory effects in Eq. (p5|) play a minor 
role only and the amplitude of motion is approximately an exponentially growing function, 
similar to the case of Newton motion from the barrier in the presence of the friction forces, 
see curve 1 in Fig. 2. The friction coefficient 7 can be derived here from Eq. (ESf) at 
ojpj i" C 1 and it is given by 7 = 7/ = Kf t = uj 2 F j 5/ r ~ r, where Upj = yjKf/Bf is the 
characteristic frequency for the eigenvibrations caused by the Fermi surface distortion effect. 
The behavior of Ag(t) is changed dramatically with an increase of the relaxation time. At 
large enough relaxation time, the descent from the barrier is accompanied by the damped 
oscillations (curve 2 in Fig. 2). These oscillations are due to the memory integral in Eq. 
(p5|). The characteristic frequency, ujr, and the corresponding damping parameter, ouj, can 
be derived from the imaginary and real parts of complex conjugated roots of Eq. ( p7|) as 
A = — ujj ± i ojr. The solution (^) takes then the form 

Ag = C c e ct + A^- Til<in sm(Et/h) + ^e"^ cos(Et/K), (32) 

where T = 2uojh and E = uo^h. In Fig. 3 we show the dependence of the instability growth 
rate parameter (, the energy of eigenvibrations E and the damping parameter T on the 
relaxation time r. 

In the rare collision regime Upj r ^> 1, the friction coefficient 7 is obtained from Eq. (p5|) 
as 7 = 7/ = Bf /t ~ 1/r. We point out that the r-dependence of the friction coefficient, 
7/ ~ 1/r, in the rare collision regime is opposite to the r-dependence of 7/ ~ r in the 
frequent collision regime. This is a consequence of the memory effects in the Fermi liquid. 
Below we will use the following extrapolation form for the friction coefficient near the fission 
barrier 



S 



7/ = Bf TTWf rf (33) 

The presence of the memory effects in the equation of motion ( p5|) changes significantly 
the trajectory of the nuclear descent from the fission barrier. The result of the solution of 

Eq. (p5|) for the trajectory Q (q) for the large amplitude motion from the saddle point qf is 
shown in Fig. 4 (solid line). The dashed line in Fig. 4 shows the trajectory obtained as a 
solution to the Newton's equation (no memory effect) 

„, , •■ dB(q) -2 dE vo Aq) 

B[q)qJr &r q = — if^~ 7/g ' (34) 

where the friction coefficient 7/ was taken from Eq. (^). In both cases we have used 



the initial conditions with q(t ) — qf, Q (t ) = y2E^o/ B j and -Ekin,o = 1 MeV and the 
relaxation time r = 4 ■ 10 _22 s. As seen from Fig. 4, the memory effect leads to the drift of 
q in time which is accompanied by the time oscillations of q along the trajectory of descent 
to the "scission" point, g sc . In Fig. 4, the time oscillations of q appear as a spiral-like 
behavior of the trajectory q (q). In both cases, the drift from the barrier is caused by the 
conservative force — dE pot (q)/dq. The oscillations appear due to the presence of the time- 
reversible elastic force in the memory integral in Eq. (f25[ ), see also Fig. 2. We point out that 
the memory integral contains the time-irreversible part also. Due to this fact, the velocity 
of the system decreases and the trajectory is shifted to the slope of the fission barrier. This 
effect is significantly stronger in the presence of the memory effects and leads to an essential 
delay of the descent process with respect to the analogous result obtained from the Newton's 
motion of Eq. fl3"4"|). The influence of the memory effect on the descent time t sc from the 
barrier to the "scission" point g sc is shown in Fig. 5. As seen from Fig. 5, in the absence 
of the memory effects (dashed lines), the descent time t sc is about 1 -f- 3 • 10~ 21 s and, as 
it should be, the value of t sc goes to the limit of non-friction motion for both the frequent 
collision regime, r —>■ 0, and the rare collision regime, r — > 00. This property of the descent 



with no-memory effects is the result of the Fermi-liquid approximation (|33"D for the friction 
coefficient 7/ in Eq. (0). In contrast to this case, the descent time t sc evaluated in the 
presence of the memory effects (solid lines) grows monotonously with the relaxation time r. 
The additional delay of the motion in the rare collision region (large r) is here caused by 
the contribution of the elastic force due to the memory integral. The elastic force leads to 



the dynamical renormalization of the adiabatic force — dE pot (q) jdq in Eq. (^) and acts 
against the force — dE pot (q)/dq. 

Let us apply our approach to the case of symmetric nuclear fission described by Eq. 
(15), assuming the Lorentz parameterization for the profile function Y(z) in Eq. (|12|) in the 
following form 

Y*(z) = (z*-e )(z 2 + e 2 )/Q, (35) 
where the multiplier Q guarantees the volume conservation, 

Q = -[Co 3 4Co 2 + C 2 2 )]/ J R 3 - (36) 
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Here all quantities of the length dimension are expressed in the Ro units. The parameter £o 
in (p5|) determines the general elongation of the figure and (2 is related to the radius of the 
neck. For ( 2 = 00 the shapes fl35|) coincide with the spheroidal ones. At finite ( 2 (C2 > for 
bound figures) the neck appears and the value C2 = corresponds to the scission point after 
which the figure is divided in the two parts for ( 2 < 0. To solve Eq. ( [15] ) we will rewrite it 
clS 8b set of two equations. Namely, 

+ e ^ m = + Mt, ?) (37) 

3=1 k=i d 9k OQi 

and 

dR^ = _R^ + j : ^ q) ^ ^ Ri{t = 0>q) = 0j (38) 
at t j=1 

where q = {91,52} = {CojC2} and the terms ~ qitjj were omitted in Eq. (|3"8"D, as the next 
order corrections. The kernel Kij(q,q) is given by 

2 r 

Kij(q,q) = -mpoVp J dr (V^V^ <\) i (r,q)) (V^V^ (r,g)). (39) 

We have performed numerical calculation for symmetric fission of the nucleus 236 U. We 
solved Eqs. (p?]) and (P8|) numerically using the deformation energy E pot (q) from Refs. |pJ,P8|. 
The scission line was derived from the condition of the instability of the nuclear shape with 
respect to the variations of the neck radius: 

f^M.O (40) 



where p nec k = C2/ \/Co(Co/5 + Cf) is the neck radius. The equations of motion (|37|) and (|38D 
were solved with the initial conditions corresponding to the saddle point deformation and the 
initial kinetic energy = 1 MeV (initial neck velocity (2 = 0). To solve the Neumann 

problem ([TJ]) for the velocity field potential we have used the method based on the theory 
of the potential, see Ref . . 

In Fig. 6 we show the dependence of the fission trajectory, i.e., the dependence of the 
neck parameter ( 2 on the elongation £ , for the fissioning nucleus 236 U for two different 
values of the relaxation time r: r = 4 ■ lCT 22 s (dashed line) and r = (dotted line). The 
scission line (dot-dashed line in Fig. 6) was obtained as a solution to Eq. (^0|). We define 
the scission point as the intersection point of the fission trajectory with the scission line. As 
can be seen from Fig. 6 the memory effect hinders slightly the neck formation and leads to 
a more elongated scission configuration. To illustrate the memory effect on the observable 
values we have evaluated the translation kinetic energy of the fission fragments at infinity, 
E'kin, and the prescission Coulomb interaction energy, -E^cmi- The value of -Ekin is the sum of 
the Coulomb interaction energy at scission point, -Ecoui, and the prescission kinetic energy 
£kin, P s- Namely, 

Eki Q =E Cou \ + -Ekin.ps- (41) 
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After scission the fission fragments were described in terms of two equal mass spheroids (see 
Ref. fl^l). We assumed that the distance between the centers of mass, d, of two spheroids is 
equal to the distance between the two halves of the fissioning nucleus at the scission point: 



d ~ 7^0^— tt! l scis • ( 42 ) 



5, Co 2 + 3CI 
4 ; °Co 2 + 5C! 

The corresponding velocity d was obtained by the differentiation of Eq. ( f42|) with respect 
to the time. The elongation, c, of both separated spheroids is defined by the condition: 

2c + d = 2Co,scis (43) 

where Co,scis is the elongation of the nucleus at the scission point. The collective parameters 
c and d and the velocity d were then used to evaluate the Coulomb energy E Co ui (see Ref. 
U) and the prescission kinetic energy E kiQjPS in Eq. ([TJ]). 

The influence of the memory effects on the fission-fragment kinetic energy, E^, and the 
prescission Coulomb interaction energy, Ecoui, is shown in Fig. 7. As seen from Fig. 7 the 
memory effects are neglected at the short relaxation time regime where the memory integral 
is transformed into the usual friction force. In the case of the Markovian motion with friction 
(dashed line), the yield of the potential energy, AE pot , at the scission point is transformed 
into both the prescission kinetic energy, E^ps, and the time irreversible dissipation energy, 
E dis , providing AE pot = i?kin,ps + -E'dis- I n contrast to this case, the non-Markovian motion 
with the memory effects (solid line) produces an additional time reversible prescission energy, 
Ef, P s, caused by the distortion of the Fermi surface. In this case, the energy balance reads 
AE pot = -Ekm.ps + -^dis + Ep tPS . We point out that the two-spheroid parametrization of the 
fissioning nucleus at the scission point given by Eqs. ( f4"2"D and (pE3"D, used in this work, leads 
to the prescission Coulomb energy Ecoui which is about 5 MeV lower (for 236 U) than the 
Coulomb interaction energy of the scission point shape Taking into account this fact 
and using the experimental value of the fission-fragment kinetic energy E^ = 168 MeV [18 



one can see from Fig. 7 that the Markovian motion with friction (dashed line) leads to the 
overestimate of the fission-fragment kinetic energy E^. In the case of the non-Markovian 
motion with the memory effects (solid line), a good agreement with the experimental data is 
obtained at the relaxation time of about r = 8 • 10~ 23 s. A small deviation of the prescission 
Coulomb energy -Ecoui obtained at the non-Markovian motion (solid line in Fig. 7) from 
the one at the Markovian motion (dashed line in Fig. 7) is caused by the corresponding 
deviation of both fission trajectories in Fig. 6. 

In Fig. 8 we illustrate the memory effect on the saddle-to-scission time t sc . In the case 
of the non-Markovian motion (solid line), the delay in the descent of the nucleus from the 
barrier grows with the relaxation time r (at r > 4 • 10~ 23 s). This is mainly due to the 
hindering action of the elastic force caused by the memory integral. The saddle-to-scission 
time increases by a factor of about 2 due to the memory effect at the "experimental" value 
of the relaxation time r = 8 • 10 _23 s which was derived from the fit of the fission-fragment 
kinetic energy E kin to the experimental value of E^ (see above). 
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IV. SUMMARY AND CONCLUSIONS 



By use of p-moments techniques, we have reduced the collisional kinetic equation to the 
equations of motion for the local values of particle density, velocity field and pressure tensor. 
The obtained equations are closed due to the restriction on the multipolarity I of the Fermi 
surface distortion up to I = 2. To apply our approach to the nuclear large amplitude motion, 
we have assumed that the nuclear liquid is incompressible and irrotational. We have derived 
the velocity field potential, 0(r,g), which depends then on the nuclear shape parameters q(t) 
due to the boundary condition on the moving nuclear surface. Finally, we have reduced 
the problem to a macroscopic equation of motion for the shape parameters q{t). Thus we 
consider a change (not necessary small) of the nuclear shape which is accompanied by a 
small quadrupole distortion of the Fermi surface. The obtained equations of motion for the 
collective variables q(t) contains the memory integral which is caused by the Fermi-surface 
distortion and depends on the relaxation time r. 

The memory effects on the nuclear collective motion disappear in two limits; of zero 
relaxation time, r — > 0, and at r — ► oo. In general case, the memory integral contains the 
contribution from both the time reversible elastic force and the dissipative friction force. 
The eigenmotion near the ground state (point A in Fig. 1) is influenced by memory effects 
through the frequency dependency of the stiffness, C'(u), and the friction, 7(0;), coefficients 
in the dispersion equation ([291). We point out that the friction coefficient in our 

approach (see Eq. Q3"I|)) changes its r-dependency from 7 ~ r in the frequent collision 
regime, Urt < 1, to 7 ~ 1/r in the rare collision regime, lurt ^> 1. Due to this fact, we 
have obtained a correct description of the zero-to-first sound transition in the nuclear Fermi- 
liquid [Q. In the limit of r — ► 00, the additional contribution (elastic force) in Eq. (25) 
appears due to the memory integral. The contribution from the elastic force is significantly 
stronger than the one caused by the adiabatic force — k Aq. The presence of the elastic force 
provides a correct A-dependence of the energy of the isoscalar giant multipole resonances. 

We have shown that the development of instability near the fission barrier (point B in 
Fig. 1) is strongly influenced by the memory effects if the relaxation time r is large enough. 
In this case, a drift of the nucleus from the barrier to the scission point is accompanied by 
characteristic shape oscillations (see Figs. 2 and 3) which depend on the parameter k of 
the memory kernel and on the relaxation time r. The shape oscillations appear due to the 
elastic force induced by the memory integral. The elastic force acts against the adiabatic 
force — dE vot (q)/dq and hinders the motion to the scission point C. In contrast to the 
case of the Markovian motion, the delay in the fission is caused here by the conservative 
elastic force and not only by the friction force. Due to this fact, the nucleus loses a part 
of the prescission kinetic energy converting it into the potential energy of the Fermi surface 
distortion instead of the time-irreversible heating of the nucleus. As mentioned above, in the 
nuclear Fermi liquid the friction coefficient 7 is a non-monotonic function of the relaxation 
time r (see Eqs. fl3lf) and (|3~3~D) providing the asymptotic behavior 7 ~ r and 7 ~ 1/r in 
both limiting cases of the frequent and rare collisions, respectively. This feature of 7 leads 
to the non-monotonic behavior of the saddle-to-scission time, t sc , as function of r in the case 
of the Markovian (no memory) motion with friction, see dashed lines in Figs. 5 and 8. In 
contrast to the Markovian motion, the memory effects provide a monotonous dependence of 
the saddle-to-scission time on the relaxation time r (see solid lines in Figs. 5 and 8). This is 
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caused by the elastic forces produced by the memory integral, which lead to the additional 
hindrance for the descent from the barrier at large r. 

The memory effects lead to the decrease of the fission-fragment kinetic energy, E kin , 
with respect to the one obtained from the Markovian motion with friction, see Fig. 7. 
This is because a significant part of the potential energy at the scission point is collected 
as the energy of the Fermi surface deformation. Note that the decrease of the fission- 
fragment kinetic energy due to the memory effects is enhanced in the rare collision regime 
(at larger relaxation time) while the effect due to friction decreases. An additional source 
for the decrease of the fission-fragment kinetic energy is caused by the shift of the scission 
configuration to that with a larger elongation parameter ( , in the case of the non-Markovian 
motion, see Fig. 6. Due to this fact, the repulsive Coulomb energy of the fission fragments 
at the scission point decreases with respect to the case of the Markovian motion. 
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FIGURES 




Fig. 1 

FIG. 1. Dependence of the potential energy E pot on the shape parameter q. 
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FIG. 2. Time variation of the shape parameter q near the saddle point B (see Fig. 1), for varios 
values of the relaxation time r. The curves 1 and 2 correspond to the values of r = 3 • 10~ 23 s 
r = 4 • 10 _22 s, respectively. 
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FIG. 3. Dependence upon relaxation time r of the characteristic energy E and width V of 
oscillations (solid lines) and the instability growth rate parameter Q (dashed line) for the curve 2 
in Fig. 2. 
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FIG. 4. Trajectory (dependence of the collective velocity dq/dt on the collective coordinate 
q) for the descent from the saddle point B (see Fig. 1). Solid line represents the result of the 
calculation in presence of the memory effects and dashed line is for the case of Markovian (no 
memory) motion with the friction forces. We have used the relaxation time r = 4 • 10 _22 s and the 
initial kinetic energy E^ in = 1 MeV. 
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FIG. 5. Dependence upon relaxation time r of the time, t sc , required to travel a nucleus from 
the saddle point B to the "scission" point C (see Fig. 1). Solid line represents the result of the 
calculation in presence of the memory effects and dashed line is for the case of Markovian (no 
memory) motion with the friction forces. The initial kinetic energy is £kin = 1 MeV. 
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FIG. 6. Trajectories of descent from the saddle point of the nucleus 236 U in the Co,C2 plane. 
Dashed line represents the result of the calculation in presence of the memory effects and dotted 
line is for the case of Markovian (no memory) motion with the friction forces. We have used the 
relaxation time r = 4 • lCT 22 s and the initial kinetic energy E^ n = 1 MeV. Dot-dashed line is the 
scission line derived from the condition (43). 
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FIG. 7. Fission-fragment kinetic energy, E^ n , (curves 1) and the Coulomb repulsive energy at 
the scission point, i?Coui> (curves 2) versus the relaxation time r for the nucleus 236 U. Solid lines 
represent the result of the calculation in presence of the memory effects and dashed lines are for 
the case of Markovian (no memory) motion with the friction forces. The initial kinetic energy is 
£ kin ,o = lMeV. 
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FIG. 8. Dependence upon relaxation time r of the saddle-to-fission time, t sc , for the descent 
from the barrier in the case of two-dimension (Co>C2) parametrization for the nucleus 236 U. Solid 
line represents the result of the calculation in presence of the memory effects and dashed line is for 
the case of Markovian (no memory) motion with the friction forces. The initial kinetic energy is 
E kin = 1 MeV. 
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